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Abstract 

Nonstationary molecular states which contain electronic coherences can be impulsively created 
and manipulated by using recently-developed ultrashort optical and X-ray pulses via photoex¬ 
citation, photoionization and Auger processes. We propose several stimulated-Raman detection 
schemes that can monitor the phase-sensitive electronic and nuclear dynamics. Three detection 
protocols of an X-ray broadband probe are compared - frequency dispersed transmission, integrated 
photon number change, and total pulse energy change. In addition each can be either linear or 
quadratic in the X-ray probe intensity. These various signals offer different gating windows into 
the molecular response which is described by correlation functions of electronic polarizabilities. 
Off-resonant and resonant signals are compared. 
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I. INTRODUCTION 


Many spectroscopic techniques involve the creation and manipulation of coherences fol¬ 
lowed by a stimulated Raman detection of a probe pulse[l-12]. The field of multidimensional 
spectroscopy was launched by looking at the delays between impulsive stimulated Raman 
events [13]. In the simplest conventional one-dimensional (ID) time-domain stimulated Ra¬ 
man technique, the molecule is first prepared in a vibrational superposition state by an 
off-resonant Raman pulse and, after a variable delay period T, the transmission change of 
a second probe pulse is detected. The transmission oscillates with T between gain and loss 
at the vibrational period, and a Fourier transform then reveals the vibrational frequencies 
[14-16]. Optical Raman techniques have been applied to study electron transfer and nona- 
diabatic dynamics at conical intersections. Using recently developed FEL and HHG sources 
[17-19], Raman techniques can be further extended to the X-ray regime [20-24] whereby the 
system is initially prepared in the superposition of valence electronic states and an X-ray 
Raman probe then reveals information about electronic, rather than vibrational, coherence. 

We had recently investigated multidimensional Raman techniques in the visible/IR regime 
[15]. Signals quadratic in the probe pulse intensity which use broadband or shaped pulses to 
probe excited state dynamics following the preparation by resonant UV actinic pulse were 
studied. Various simulation protocols which differ by the level of complexity and computa¬ 
tion cost were compared. Starting with the most computationally expensive expensive exact 
numerical propagation, we developed more approximate sum-over-states expansion and the 
Stochastic Liouville Equation (SEE) to describe complex dynamics. Excited state prepa¬ 
ration has been limited to excited electronic state populations which do not carry phase 
information. In our recent work [25] we had extended the SEE approach to both linear and 
quadratic off-resonant Raman signals with shaped pulses. The Raman resonance line shapes 
were calculated for systems with nonequilibrium initial state prepared by a resonant actinic 
pulse. We then investigated the relative phase dependence of time evolving Stokes and anti- 
Stokes components of the spectra in quadratic signal with shaped pulse [26]. The results 
show the energy redistribution between different spectral components in Raman signals. 

This paper presents several electronic Raman techniques performed with X-ray pulses that 
carry higher levels of information [21, 27, 28]. The same techniques apply to conventional 
vibrational Raman as well. Using an intuitive diagrammatic approach, we derive expressions 
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for three basic detection protocols of the probe: (i) the change in number of photons 
(ii) the change in the transmitted probe energy and (Hi) the frequency-dispersed probe 
transmission In addition, signals may be linear or quadratic in the probe intensity 

and may utilize broadband or hybrid shaped (a combination of broad and narrow) probes 
as well as resonant or off-resonant with material transition. The energy exchange between 
held and matter is discussed. In the case of an off-resonant probe, the number of photons 
is conserved and = 0. For a hybrid pulse composed of a narrowband and a broadband 
component, we see oscillations between Stokes/anti-Stokes components [26] with the delay 
T, whereas for a single broadband pulse the entire pulse envelope oscillates, smearing out 
the spectral features in the hybrid signal [29]. 

We present a systematic organization scheme for stimulated Raman multidimensional 
spectroscopies applicable when the preparation of some nonstationary state, the nature or 
dynamics of which is the object of study, is temporally well-separated from the detection 
process. Many conventional spectroscopic techniques (e.g. pump-probe type) fulhll these 
conditions [30-32]. The scheme takes the inital state as a given and classihes signals by 
their dependence on detection parameters. These parameters control which features of the 
nonstationary state are observable and how they manifest. In particular, we examine dif¬ 
ferent choices for the held (whether the held spectrum is resonant or oh-resonant with the 
material transitions), the intensity scaling with the detecting held (linear or quadratic), the 
spectral shape of the detecting held (broadband or hybrid broad-narrowband). These are 
only a few of the possible combinations of parameters that could describe spectroscopic 
detection of a nonstationary state. Similar ideas were pursued for oh-resonant X-ray scat¬ 
tering (dihraction) [33, 34] and for spontaneous emission following impulsive Raman X-ray 
excitation [35]. 

While oh-resonant Raman signals are simpler to analyze, resonant pulses are more se¬ 
lective to a given atomic core transition and provide additional specihc information about 
molecules. We examine diherent simulation protocols and derive expressions that may be 
used for direct numerical propagation of the wavefunction in Hilbert space, which includes 
all degrees of freedom explicitly, or the density matrix in Liouville space [36], which gives a 
simple picture of bath ehects and allows for a reduced description via the Stochastic Liouville 
Equations (SLE) [15, 25, 37]. 

In the oh-resonant regime the held-matter interaction Hamiltonian is the product of held 
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intensity and a molecular polarizibility, while in the resonant regime it is the dot product of 
held amplitude and molecular transition dipole. The off-resonant signals are explored first 
and expressions for each of the three detection protocols for signals linear and quadratic in 
the field intensity are derived. We then present plots of these signals for a simple model 
system and discuss the properties of each signal and what information it reveals. This 
analysis is then repeated for the resonant signals using the same model. We conclude with 
a comparison of off-resonant versus resonant signals and a general discussion of the utility 
of this signal classification scheme. 


II. DETECTION PROTOCOLS FOR OFF-RESONANT RAMAN SIGNALS 

-X 

- y 


- d 

c 

- a 

FIG. 1: Level scheme for the broadband X-ray Raman experiment, a, c, d are 
low-frequency, valence electronic excitations and x, y are high-frequency core excitations. 

We will investigate stimulated X-ray Raman signals from a molecular model system with 
a band of valence states a, c, d, etc. and a core excited state band x, y as shown in Fig. 1 
(though, as mentioned in the introduction, the formalism applies as well for other different 
energy regimes e.g. x, y valence excitations and a, c, d vibrations). We assume that the 
system is prepared at time tq in a superposition of valence states ipi and is monitored by 
interaction with a broadband probe pulse centered around time to- This is depicted diagra- 
matically in Fig 2 for both linear and quadratic signals. Note that, although off-resonant 
diagrams can be deduced from the corresponding resonant diagrams, they represent pros- 
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FIG. 2: (Color online) Loop diagrams representing the off-resonant linear (a) and 
quadratic (b) signals as well as the resonant linear (c) and quadratic (d) signals. In all 
diagrams, the system is assumed to be prepared from the ground state by some unspecihed 
process (depicted by the grey rectangles). The preparation process terminates at time tq 
after which the system evolves freely until it begins interacting with the probe which is 
taken to be centered at time to- The delay parameter T = to — tq is therefore shown next 
to the diagrams. Note that in the off-resonant case, no core-valence coherences are created 
while these are created for certain times in the resonant case. For the linear processes ((a) 
and (c)), Uca > 0 implies a red (Stokes) contribution while the reverse condition implies a 
blue (anti-Stokes) contribution. For quadratic signals, this analysis holds only for one of 
the two relevant diagrams and so the quadratic signals are more difficult to interpret (since 
it depends on the state d). Details of the loop diagrams are given in [38]. 
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esses governed by different interaction Hamiltonians (Eq. (1) versus Eq. (26) for resonant 
interaction). We therefore include both sets of diagrams for completeness. A brief summary 
of different detection protocols for linear and quadratic signals and their most striking fea¬ 
tures in the off-resonant case are shown in Table I. The held-matter interaction Hamiltonian 
for an off-resonant Raman process in the interaction picture is given by 

//'(«) = «'">(*) If Wt. (1) 

where represents the off-resonant polarizability, which is a Hermitian operator in the 
valence space = Q;f-|-Q; is real). In the following, we assume that the electric held consists 
of short pulses which are temporally well-separated from the preparation process and Fourier 
transformable S(t) = This simplihes the analysis while retaining all the 

essential physics. The frequency-dispersed transmission of the probe pulse (heterodyne- 
detected, frequency-dispersed photon-number change) is given by 

9 ^OO 

= -X j dte*"(*-*°)r(a;)T(t - to)((/|ai°H^)e-^ (2) 

where pi = \'ipi)(ppi\ is the initial density matrix (immediately following the end of the state 
preparation process), double brackets denote Hilbert space operators that are vectors in 
Liouville space, ((/| represents the trace operator in Liouville space, and X{1Z) denotes 
the imaginary (real) part. We adopt a superoperator notation that provides a convenient 
bookkeeping of time-ordered Green’s functions. With any ordinary operator A we associate 
two superoperators dehned by their action on an ordinary operator X as Al = AX (action 
from the left) and Ar = XA (action from the right). We further dehne the symmetric and 
antisymmetric combinations A_|_ = \{Ar + Aj?), A_ = \{Al — Ar). The H'_ exponential 
represents evolution of the matter with the probe held (since the preparatory processes 
have ceased) until the hnal interaction with the probe. We can then expand the signal 
perturbatively in H'_ to obtain signals with the desired held-scaling (linear or quadratic in 
intensity are considered in this paper). 

Alternatively one can detect the total change in the photon number which is given by 
the zeroth spectral moment of 

S''^\to,To) = f (3) 

J-oo 27r 
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Summary of Off-Resonant Stimulated Raman Techniques. 

Signal 


S(E) 

Slb 

Eq. (10): Oscillatory gain/loss pattern 

shows Stokes/anti-Stokes oscillations in 

T. 

Eq. (12): No spectral resolution due 

to u) integration. May visualize weak 

transitions due to weighting ujac- 

Slh 

Eq. (13): High spectral resolution com¬ 
pare to Slb- Each peak oscillates at it’s 

frequency ojac and with phase (pac- 

Eq. (15): Does not carry new informa¬ 
tion compared to S^d<i\ 

to 

Eq. (22): Always Stokes due to domi¬ 
nance of populations (which do not con¬ 
tribute to linear signals). 

Eq. (23): Transition spectra can be 

achieved by Eourier transform over T. 

Sqh 

Eq. (24): spectral and temporal resolu¬ 
tions are higher than Slh (not conju¬ 
gated) but retrieval of frequencies and 

phases is more complicated. 

Eq. (25): with Eourier transform over 

T can visualize weak transitions. 


TABLE I: Summary of off-resonant Stimulated Raman Techniques. Note, that signals 
vanish for both linear and quadratic detection protocols, so it has been omitted. 

This gives 

9 ^OO 

S^^\to,To) = -I j dt\S{t-to)\‘^{{I\af\t)e--^^^-^^^'^^\pi{To))). (4) 

A third type of signal is given by the total change in the energy of the transmitted pulse 
which is the hrst spectral moment of 

S^^\to,To) = j ^lTwS^^‘^\u,to,To), (5) 

also given by 

/ OO 

dtS* {t - to) S {t-to) {{I \ (^) e" ^ I Pi (to) )). (6) 

•OO 

Sif<^) is the most detailed measurement whereas and are its zeroth and hrst spectral 
moments. The latter two techniques are simpler to implement experimentally since they do 
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not require frequency-selective detectors or spectrometers. In the following we calculate the 
three signals (2) - (6) for a model system probed by different choices of pulse shapes and 
parameters. 


A. Off-Resonant Signals Linear in the Probe 

This technique is essentially a stimulated Raman pump-probe and is the off-resonant 
analogue of transient absorption. Here, the photons undergo a redistribution among held 
modes rather than being absorbed. Expanding Eq (2) to linear order in the held, we obtain 

/ OO 

- To)\p.)), (7) 

■OO 

where Q{t) = is the Liouville space Green’s function of the molecule. In 

some applications, it is necessary to carry out an exact numerical propagation of the wave 
function rather than using the density matrix. We can then set 

{{I\a^L^Q(t - To)\pi)) = - ro)|V’i), (8) 

where G{t) = {—i/h)9{t)e~^^°^ is a Hilbert space Green’s function. This result can be 
further recast in the Heisenberg picture as 

- To)\pi)) = - To)\'ipi). (9) 

1. Linear Broadband (LB) Probe 

We now expand the signals (9), (3), and (5) in molecular eigenstates . Using Fig. 2(a), 
we obtain 

Slb\^^T) = -^^|T(a;)||T(a;-a;„c)|a(°Vac|sin0„e(T), (10) 

a,c 

dfl(r) = o, (11) 

a,c 

8 


( 12 ) 



Here, pac = |pac|e*'^“‘= is the initial density matrix, the phase is given by 4>ac{T) = UacT — 
and we neglect the linewidth for the valence states. In Eq. (10), terms in which ojac > 0 (cuac < 
0) lead to broad peaks below (above) the central pulse frequency. We denote these as “red” 
and “blue” peaks respectively. Each pair of states therefore generates a complementary pair 
of red and blue contributions. The two oscillate in T with a tt phase shift (as sine is an odd 
function of its argument). The frequency-dispersed signal therefore oscillates between Stokes 
(positive red contributions and negative blue) and anti-Stokes (negative red contributions 
and positive blue) processes as shown in Fig. 3(al). 

Equation (12) is plotted in Fig. 3(bl) but, aside from indicating that some nontriv¬ 
ial phase is involved (since the signal doesn’t vanish for T —)■ 0), not much information 
is directly apparent from the time-domain measurement. Taking the Fourier transform 
f dTe^^'^S{T) = S'(r2) (Fig. 3(cl)) reveals peaks at each transition energy Uac whose heights 
are given by the factor ujac(^ca\pac\- The ratios of these peak heights thus give information 
on the polarizabilities or the magnitudes of the initial coherences. Finally, it is worth noting 
that (paaiT) = 0,VT so that populations do not contribute to the off-resonant linear broad¬ 
band signal. These signals therefore provide a baekground-free deteetion of the electronie 
coherenees. 

2. Linear hybrid (LH) probe 

A hybrid probe is a shaped pulse consisting of a broadband attosecond pulse S{u) = 
|To(ci;)|e*'^“ and a narrowband femtosecond pulse Si(t) = centered at frequency 

Ui- In this case, Eq. (7) yields 

= -y|To(a;)||Ti|^a(°Vac|sin0i^(T)(5(a;-a;i -Wac), (13) 

a,c 

where p\cO^) = oJacT — + po — Pi- Note that if the broadband and the narrowband 

components have the same phase, po = pi, then Pl^iT) = pac(T). Due to the dependence 
on the relative phase of the pulses po — pi, observation of the linear hybrid signals requires 
phase-control (averaging over random pQ, pi causes the signal to vanish). Eq. (13) yields 
clearly resolved Raman resonances unlike the LB case where the red {uac > 0) and blue 
{uac < 0) components only enter through the pulse envelope. This signal is depicted for 
oji = ojq in Fig. 4(al), shows sharp peaks at each transition frequency Uac- Of particular 
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significance is the ability to extract the phase 0^^ from the oscillations of the separate peaks. 
This information is not available in the LB signal. 

The photon number signal is slightly more complex: because the hybrid pulse contains 
both broad and narrowband components we need to take into account the contribution to 
the signal where the last interaction is with narrowband component. The total change in 
photon number vanishes, as in the broadband case, 

S'">(i..„T)=0. (14) 

The transmitted energy change of the shaped pulse similarly contains both narrowband and 
broadband components and the total energy change of the shaped pulse is given by 

= -y ^|Ti||To(a;i - a;ac)|ttiaVac|ha;acSin0^^(T). (15) 

a,c 

Note the similarity between this and Eq. (10) (they only differ by the factor ujac£{,^))- The u 
integration erodes the sharp resolution afforded by the narrowband pulse, leaving virtually 
the same result as the signal (Eq. 10). The signal is shown in Fig. 4(bl). Note 
that the presence of the ojac factor which reverses the sign under a ^ c relative to Eq. (10) 
rendering both red and blue contributions the same sign for a given T. The sign of 
then indicates whether the process is Stokes < 0) or anti-Stokes (S'^^^(T) > 0). 

Energy conservation implies that the pulse energy change and the molecular energy change 
must have equal magnitude and opposite sign. In particular, the diagram in Fig. 2(a) 
assumes that, after preparation, the molecule is in the superposition of states a and c and 
the hnal state of the molecule after interacting with transmitted pulse is a population cc. 
Similarly using the permutation a <H- c one can end up in the final state aa. If the energy of 
state a is higher (lower) than that of c this is a Stokes (anti-Stokes) process. The molecular 
energy change is given by 

S^^\to,To) = - '^hWac{paa{to,ro) - pcc{to,ro)). (16) 

a,c 

This result holds for an arbitrarily prepared molecular state prior to the probe. 

B. Off-Resonant Signals Quadratic in the Probe 

The quadratic signal (i.e., second-order in probe intensity) is obtained by expanding the 
exponent in Eqs. (2)-(5) to hrst-order. The signal can be read off the two diagrams shown 
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in Fig. 2(b), yielding 


r*oo pt 




' —OO k/ —OQ 

S^)n(+ _ 


X ((^l«L - t')a\’Q{t' - ro)|pi)), 


( 17 ) 


S™(<o.To) = 12! / dt j dt'\S(t - ta)\‘^\S(t' 

J —OO J —OO 

X {{I\aPg{t - - To)\pi)), 


(18) 


S^Q^\to,To) = 2hI 


7 *CXD pt 

dt / dt'£*{t— tQ)£{t— to)\£{t' — to)Y 


' —OO J —OO 


X {(/|a5°y(i - - To)|pi>), 


(19) 


The corresponding Hilbert space expressions which are suitable for wavefunction-based sim¬ 
ulations are 


{{I\aY’g(t - t')a^yg(t' - To)\pi)) = -ih[(^/’i|G^(t - ro)a^°^G(t - - ro)|^/'i) 

-{ipi\GHt' - To)a'''lGl(! - - To)!!/!!)] 


( 20 ) 

Eq. (20) may be alternatively recast in the Heisenberg representation 

{{I\af'^g{t - t')a^°^g{t' - To)\pi)) = -ih[{gi\a^^\t - To)a^^\t' - To)\gi) 

- To)a^^\t - To)\gi)]. (21) 


1. Quadratic Broadband (QB) Probe 


The frequency-dispersed quadratic signal (17) for a broadband probe expanded in eigen¬ 
states reads 

dui I 




'QB 


■ I^ 11^ (^1) I «cd «ia IPac I COS (T) 

a,c,d 


J 271 

X ||f(ip - x>*)|i'(x>i +x>i„)| - |£(x) + x)j„)||£(x)i - x>*)||, 


( 22 ) 


where the first (second) term in the square brackets corresponds to the left (right) diagrams 
in Fig. 2(b). Eq. (22) has a more complex dependence on the Raman shift Uac than the 


11 



linear signal (10). The field envelopes are now shifted by the electronic transition frequencies 
ojad and Udc which yields difference Udc — ^da = ^ac- The quadratic signal oscillates with a 
phase that depends on states d other than a and c that create the resonance and involves 
the phases of the polarizability 0 ;*'°^ 

Comparing this to the corresponding linear signal (Eq. (10)) shows that a 7r/2 phase 
change. As a result, the fact that 4>aa{T) = 0 no longer eliminates contributions from 
populations and they form a time-dependent background to all quadratic signals. Because 
the initial density matrix is assumed to be perturbative (so that the ground state dominates 
the populations) the contribution due to populations is primarily Stokes overall. Since the 
contribution due to the populations is Stokes and the oscillating coherences are too weak to 
overcome this, the overall process is Stokes at all times T (as illustrated in Fig. 3)(a2). 

The integrated photon number vanishes = 0. The energy change signal (5) is 

given by 

4B(r) = -4/^/+ 

a,c,d 

X {hWda + f^dc)\Pac\a^S<=^d(! COS0ac(T). (23) 

Since the populations contribute a static off-set, there is now a strong zero-frequency peak 
in the Fourier transform of Eq. (23). Additionally, the peak heights are no longer as sim¬ 
ply related to the aac, Pac as in the linear case. In particular, since the peak heights are 
determined by a free summation over the intermediate state d, the terms in this sum can 
interfere constructively or destructively leading to enhanced or suppressed peaks (note that 
the Uac = 4 eV peak is suppressed in Fig. 3)(c2). 

2. Quadratic Hybrid (QH) Probe 

By expanding the quadratic frequency-dispersed transmission signal (17) in eigenstates 
we obtain 

= -^|To(x;)||Ti|2^4°jaSVac|cos0„c(T) 

a,c,d 

X [|To(c(;i -|- u:da)\d{io — Ui — Udc) — |To(ci;i — u:dc)\5{u — u}da)] ■ (24) 

As in the linear case, the narrowband pulse allows us to clearly resolve the transition peaks. 
Unlike the linear signal (13), the quadratic hybrid signal is independent of the phases of 
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the narrowband and broadband pulses 01 and 0 o and is therefore observable without phase- 
controlled pulses. Just as in the broadband case, the signal is sensitive to populations which 
contribute a static Stokes spectrum. In Fig. 4(a2), we plot this signal for ui = ojq as well as 
separate plots for the contributions due to populations and coherences. 

To calculate the integrated signals we must also include contributions from diagrams 
whereby the last interaction is with narrowband pulse £i. The total photon number change 
then vanishes = 0. The corresponding pulse energy change which includes both 

broadband and narrowband components is 

Sqh{^ 1,T) = \pac\ cos (paciT) 

a^c,d 

X [hWdc\£o{<^l + UJdc)\\£Q{(^l + i^da) \ + huJda\£o{(^l “ | “ i^dc)W- (25) 

Note that in the quadratic case the energy flux involves also another state d in addition to 
states a and c, so the different fluxes corresponding to all relevant pairs of states should be 
added to get the overall energy change of the pulse. 

C. Discussion of Off-Resonant Signals 

In the previous section we considered a variety of heterodyne-detected signals in the off- 
resonant regime. Here, we present these signals calculated for a simple model consisting of 
a ground state and two valence excitations at ^^(c) £ {0,2.7, 6 . 8 } (eV). The polarizabilities 
aac are symmetric but otherwise random numbers of 0(.l — 1 ) and the pac are random but 
taken from a pure perturbative state (i.e., a state 10 ) in which the ground state amplitude 
is near unity in magnitude). 

Figure 3 shows the linear and quadratic off-resonant broadband signals side-by-side. We 
do not plot since it vanishes in the off-resonant regime. For the frequency-dispersed 
signal (Eqs. (10), (13), (22), and (24)), the broadband detection makes discerning individual 
transitions impossible but nicely illustrates the Stokes/anti-Stokes oscillations in time. In 
the linear case, the red contributions are positive and the blue are negative at T = 0 so 
the process is Stokes at this time and oscillates with T. The quadratic detection is always 
Stokes due to the effect of populations. The time-domain energy signal (Eqs. (12), (15), 
(23), and (25)) reveals again that the quadratic process is always Stokes while the linear 
process oscillates between Stokes and anti-Stokes. Transforming these signals reveals the 
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FIG. 3: (Color online) (al) and (a2): Off-resonant linear and quadratic 
frequency-dispersed signals (Eqs. (10) and (22)) at various times T after state preparation 
(T advances in units of ~60 attoseconds as the dashes lengthen), (bl) and (b2): 
Off-resonant linear and quadratic broadband energy signal as a function delay time T 
(Eqs. (12) and (23)). Note that the linear signal oscillates about zero but the quadratic 
signal has a static offset corresponding to the contribution from populations, (cl) and (c2): 
Fourier transforms of (bl) and (b2). In the linear case, peaks corresponding to all cuca 
coherences are visible. In the quadratic case, there is a large central peak at zero 
corresponding to the populations. In the quadratic signal, the states c — a are coupled 
indirectly through a third state d and the summation over intermediate states can suppress 

or enhance the Raman peak magnitudes relative to their linear proportions. 
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FIG. 4: (Color online) (al) and (a2):Off-resonant linear and quadratic hybrid 
frequency-dispersed signals (Eqs. (13) and (24)) at various times T after state preparation 
(T advances in units of ~60 attoseconds as one goes up the vertical axis). The quadratic 
signal is split into the static contribution due to populations (a2 top spectrum) and the 
time-dependent contribution due to coherences (a2 lower spectra). In contrast to the 
broadband signals, the narrowband pulse allows us to resolve the individual transition 
peaks, (bl) and (b2): Off-resonant linear and quadratic hybrid energy signals (Eqs. (15) 
and (25)) as a function of narrowband frequency uji and delay time T (T advances in units 
of ~60 attoseconds as the dashes lengthen). The broadband detection renders individual 
transition peaks unobservable and the result is similar to the broadband case (Eq. (23)) 
but symmetric (rather than anti-symmetric) about u = cjq. Note that Sqj^(T) < 0 
revealing that the quadratic process is always Stokes while the linear process oscillates 

between Stokes and anti-Stokes. 
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ojca transitions (though the quadratic possesses a large hi = 0 peak due to populations that 
is of course missing from the linear signal). 

For comparison, the linear and quadratic hybrid signals are shown in hgure 4. The 
spectral resolution of the narrowband pulse gives sharp peaks at the u — ui = Uac- In the 
linear case, each Uac peak is due to a term that oscillates with this same frequency. However, 
in the quadratic case, each oscillating term contains peaks at all different u^a, ^dc (as per 
Eq. (24)) and the phases of the peak oscillations do not as directly reveal the phases 0^^. 
Note that the free summation in the quadratic signal implies that peaks will be visible even 
for states not initially occupied (i.e. for which pac = 0) while the peaks observable in the 
linear signal are restricted to those initially occupied. 

In summary, shows an oscillatory pattern of gain and loss features in the red and 
blue regimes that depend on the initial phase of pac- This information is integrated out 
in where the entire probe pulse envelope exhibits Stokes (loss) or anti-Stokes (gain) 
shifts periodically [39]. The broadband signals alone do not possess a sufficient resolution 
to directly observe the transition spectra and only allow us to access it through the Fourier 
transforms of the energy signals Utilizing a hybrid broad-narrow pulse combined 

with frequency-dispersed detection allows high spectral and temporal resolution thus permit¬ 
ting spectral snapshots to be taken that clearly resolve all transitions and therefore permit 
extraction of the phases from the oscillation patterns. In the linear case, this requires control 
of the relative pulse phases (the phases cancel in the quadratic case). The quadratic and 
linear signals generally carry the same information about coherences but this is accessed 
without the background due to populations in the case of linear signals. One important 
caveat to this is that the magnitudes of the various transition peaks can be enhanced or 
suppressed by the interference of different pathways (characterized by the intermediate state 
d). Since each oscillation term contributes to multiple peaks through this intermediate state, 
a given peak Uac will not oscillate at the corresponding phase 0^^, making retrieval of the 
wavepacket’s phase more difficult in the quadratic case. 

Equations (9) and (20) can be used for the numerical simulation of signals from systems 
undergoing complex dynamics. The Liouville space form can describe the reduced density 
matrix by using the SEE. The Hilbert space form is suitable for direct numerical simulation 
of the wavefunction that includes all degrees of freedom. It can apply e.g. to nonadiabatic 
dynamics at conical intersections. The eigenstate expansions (Eqs. (10)-(15) and (22)-(25)) 
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may be used when the probe is impulsive so that the eigenstates do not vary during the 
probing process. Appendix A extends this treatment of off-resonant signals to a two-pulse 
time-domain experiment and considers the preparation by an actinic pulse. 

III. RESONANT STIMULATED RAMAN TECHNIQUES 

The off-resonant techniques considered in section II conserve the number of probe pho¬ 
tons. Photon energy gets redistribnted among the broadband modes: vanishes while 

SiE) 

is hnite. When the X-ray pulses are resonant with core transitions, true photon ab¬ 
sorption can take place. This renders finite as well revealing new matter information. 
Below we discuss both linear and qnadratic resonant signals. 


A. Linear Probe 

The signal linear in the probe intensity is given by the diagram shown in Fig. 2(c). Dne 
to resonant excitation, the held-matter interaction Hamiltonian may no longer be recast 
using a freqnency-independent, off-resonant polarizability bnt should rather be described by 
a dipole interaction Hamiltonian in the rotating wave approximation (RWA) 

= V{t)E^{t) + V^{t)E{t). (26) 

The freqnency-dispersed signal can be read off the diagram 

/ OO pt 

/ dt'£{t' -to){{i\VLg{t --To)\p^)). 

■OQ J —OO 

(27) 

The integrated photon number is given by 

/ OO pt 

dt / dt'8*{t - - to){{I\VLgit - t')vlg{t’ - To)\p^)). (28) 

■OO J —OO 

The pulse energy change is given by 

s[^\t,,To) = 2h^n r dt f ^°^ T(F - to)((/|f4g(t - t')vlg{t' - To)\pg). 

J— OO J— OO (liQ 

(29) 
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The signals (27) - (29) can be alternatively recast in the Hilbert space form suitable for 
numerical propagation by setting 

mVLQ{t - - ro)|A)) = - To)VG{t - t')V^G{t' - toM) (30) 

These signals are often referred to as “transient absorption”. 


1. Linear Broadband (LB) probe 


Expanding the frequency-dispersed transmission (27) as a sum over states and evaluating 
time integrals we obtain 

Slb\^:T) = \pac\ (^aif^^'(a;) sin0ac(T) - cos0ac(T)) , (31) 

a,c 

where ^ac(T) = UacT — 0^^) a” represent the real and imaginary parts of the polar¬ 

izability 




E 


(^)^fc(^ T T ^ca){L^k ■ ■ Vxc) 

UJ -\- UJj — Uxa + ^Ta; 


(32) 


where Um, m = j,k is the central frequency of the pulse given by Sm{t) = JY 
where we redehned the frequency domain amplitude ^ 0 ,( 0 ;) to be centered at zero frequency 
(for a Gaussian pulse with bandwidth am, Tm(ca) = ^ ^), Tg is the inverse excited 

state lifetime, which is assumed to be shorter than states a and c. Thus, the polarizabil¬ 
ity implicitly depends on the pulse parameters, such as its central frequency. For j = k, 
Eq. (32) reduces to the commonly used polarizability (see Eq. (5) of [27]). Note that in the 
off-resonant case the polarizability matrix elements are real and the second term in Eq. (31) 
vanishes and the Stokes and anti-Stokes components oscillate with opposite phase as seen 
in Section II. The integrated photon number (28) is 


dslr) = ^ E COS'#>ac(r), 


(33) 


a,c 


where the integrated polarizability is given by 




(34) 


and we used the symmetry ~ signal is given by the 

imaginary part of the polarizability. Therefore, when the polarizability is real (as in the 
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off-resonant case) the signal vanishes, which is consistent with our earlier result of Eq. (11). 
Again the Stokes {uac > 0) and anti-Stokes {uac < 0) components oscillate with an opposite 
phase. The energy change of the pulse signal (29) is given by 




sin- (/3£*’" + /3if'")cos.#'„c(r) 


(35) 


a,c 

where the tensor represents the hrst moment of the polarizability: 

/3a''> = /_“^fi(n> + n.,)y?>(in). (36) 

Unlike the photon number signal (33), the energy change (35) involves both the real part of 
the polarizability tensor and the imaginary part of its first moment. Eq. (35) reduces to Eq. 
(12) in the off-resonant case, where the polarizability is real and the second term vanishes. 


2. Linear Hybrid (LH) probe 

The frequency-dispersed transmission of the broadband pulse £{uj) (27) in the presence 
of the narrowband pulse £i (hybrid probe) is given by 

Slh\^ = -^^\pac\ (u) sin (j)ac{T) - a),°^^"(a;) cos0ac(T)) , (37) 

a,c 

where the ui dependence is now implicitly included in by setting Si{u) = SiS{u). The 
integrated photon number (28) for a shaped pulse has to be calculated differently than the 
broadband case. In particular, one has to take into account the photon number change in 
both the broadband and narrowband pulses. This yields 

sin0ac(T) - cos0ac(T)) , (38) 

a,c 

where the hrst (second) terms in each square bracket represents the last interaction with 
the broadband (narrowband) pulse. For the hybrid pulse, both a' and a” contribute to the 
signal, thus providing additional molecular information than the simple broadband pulse. 
Furthermore, the form of the anti-Stokes and Stokes polarizabilities suggests that the signal 
depends exclusively on the phase difference between both broadband and narrowband helds 
00 ~ 01 which provides an additional control knob. Note, that in general the hybrid polar¬ 
izabilities are not symmetric under permutation of their indices. However in the limit 
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( 39 ) 


when 00 = 01, we have a'^y^ = afx'* and Eq. (38) yields 

S^Lh{^1^T) = ldac|[aia^^" + COS (j)ac{T). 

a,c 

Finally, the total pulse energy change (narrowband and broadband components) can be 
obtained from Eqs. (38) 

= sin0ac(T) - + 0£°)"] COS0ac(T)) . (40) 

a,c 

B. Quadratic Probe 

The diagrams for the quadratic signal are depicted in Fig. 2(d). We read the signal 
from the diagrams in Hilbert space (corresponding Liouville space expressions can be easily 
derived by expanding the loop diagram in the set of ladder diagrams) and obtain 

/ OO pt POO pt\ 

dt dt' dti dt[S*{u)S{t' 

OO J—OO J—OO J—OO 

X [s*{h - - to)(0*|Gt(t - ro)HG'(t - t’)V^G{t' - h)VG{h - t'^)V^G{t\ - ro)|00 

+ (T(H - to)£*{t\ - - to)VG\g - t\)V^G\t - G)VG{t - t')V^G{t' - ro)|00], 

(41) 

where the hrst (second) term in Eq. (41) represents the left (right) diagram in Fig. 2(d). 
The corresponding photon number and energy change signals are given by Eqs. (3) and (5) 
respectively. 


1. Quadratic Broadband (QB probe) 


For a broadband pulse we expand the signal (41) in eigenstates which yields a compact 
formula 


Sgsi^^T) = -x|j^[a2^^(a;)4f^ - ^(x;)4f ^>ac(T), 




(42) 


a.c.d 


where Pac{T) = \pac\e in the off-resonant case a" = 0 and Eq. (42) reduces to 

Eq. (22). The integrated photon number signal (3) reads 


= -^Yl sin0a,(T) - afG" cos0ac(T)], 


(ssyrissy 


(ssy 




(43) 


a,c,d 
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which clearly vanishes in the off-resonant case when a” = 0. The total energy change of the 
pulse (5) is given by 




a,c,d 




(44) 


2. Quadratic Hybid (QH) probe 


The frequency-dispersed transmission of the broadband component of a hybrid broad- 
narrow probe (Eq. (41)) reads 




.2i 


(oi)o o.,(io) .,(01)/ o.,(oi)q 




(45) 


a,c,d 


which is similar to Eqs. (9)-(10) of [28] for the model with zero linewidth of electronic states. 
The integrated photon number signal (3) which includes the change in both broadband and 
narrowband helds is 

c(^)^ (01) (10) I (10) (01) (01) (01)* (10) (io)*i 

T) = -X— ]pac{T) (46) 


a,c,d 


and the energy change of the pulse (Eq. (5)) is given by 

, Ti = -T— ^ - /q( 01)^(01)* _ /Q(10)^(10)*n /^x /x^x 

?j4 / ^ \-Pcd ^da ' Pcd ^da Pda ^dc Pda ^dc JPacl^ )■ 


a,c,d 


C. Discussion of Resonant Signals 


For the simulations of the resonant signals, we extended the model used in the off- 
resonant case to include core states with energies ex{y) € {136,141.5,149.5} (eV). Figure 5 
shows the resonant linear and quadratic broadband signals. The frequency-dispersed signals 
(a/b) contain the valence-core resonances but not the Raman (this is due to the lack of 
held resolution). However, the Raman transitions Uac are in the same region as the shifted 
valence-core transitions u^a — so they appear similar. Since the u^a peaks can arise from 
any a, c pair in the summation, the phases 0^^ are not directly accessible as the phases of 
peak oscillations. Taking Fourier transform of the energy of photon number signals (c) and 
(d) gives peaks corresponding to the Raman transitions only since these are the oscillation 
frequencies. As in the off-resonant case, the linear and quadratic signals have different 
relative peak intensities controlled by the different forms of the coupling. 
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Figure 6 shows 2D spectra of the resonant linear (left column) and quadratic (right col¬ 
umn) hybrid frequency-dispersed signals for populations only (top panels) and populations 
and coherences together at different times (remaining panels). As Ui — cuq varies, peaks 
corresponding to Raman transitions (uJca) move along diagonals forming resonant streaks. 
In contrast, the core transitions do not vary with tui and therefore, each transition forms 
a series of repeated vertical peaks where the cj^a transition intersects the diagonal streaks. 
The magnitude of these peaks then reveals the strength of the coupling between the core 
and Raman transitions. Note that the population contribution makes only a single diagonal 
streak in the linear case but results in all streaks in the quadratic case (this is due to the 
summation over the intermediate state d). For the linear signal, we can therefore identify 
each linear streak with a particular a, c pair in the summation and each streak will un¬ 
dergo Stokes/anti-Stokes oscillations at the respective phase <hac(F), therefore allowing the 
retrieval of the phases 0^^. 

Figure 7 displays linear and quadratic hybrid energy and photon number signals. In both 
cases, the contribution due to populations is stronger than that from coherences. This is 
due to the resonant nature of the signal since the population contribution vanishes for linear 
off-resonant signals. Unlike in the off-resonant case, the signal contains sharp peaks due to 
the valence-core transitions but the same lack of held resolution prevents direct identihcation 
of the Raman transitions. These can however be obtained via the Fourier transform. Figure 
8 shows the Fourier transform (magnitude) with respect to delay time T of the linear and 
quadratic photon number and energy signals. The quadratic is notably weaker but all four 
signals show the same basic pattern of uj^a — resonances along the y-axis and D = ujac 
resonances along the x-axis. In Appendix C, we consider the resonant signals in the context 
of a nonstationary state prepared by an actinic pulse. 


IV. CONCLUSIONS 

We presented a systematic classihcation scheme applicable to spectroscopic techniques 
in which a nonstationary state’s creation is well-seperated from it’s detection (pump-probe 
style spectroscopies). In this scheme, a particular spectroscopic technique is specihed by 
the choice of state preparation and detection procedures. We consider various detection 
procedures and examine what sort of information one may obtain and in what ways the 
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FIG. 5: (Color online) (a): Resonant linear broadband frequency-dispersed signal for the 
populations (time-independent) and coherences at T = 0 (left) and time evolution of 
coherences (right), (b): Resonant quadratic broadband frequency-dispersed signal for the 
populations (time-independent) and coherences at T = 0 (left) and time evolution of 
coherences (right). For coherences in both (a) and (b), T advances in units of ~60 
attoseconds as one goes up the vertical axis, (cl) and (c2): Resonant linear (cl) and 
quadratic (c2) broadband energy signal !®^rier transformed, (dl) and (d2): Resonant 
linear (dl) and quadratic (d2) broadband photon number signal Fourier transformed. 
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FIG. 6: (Color online) Resonant linear (left column) and quadratic (right column) hybrid 
frequency-dispersed signals. The top of each column is the time-independent contribution 
due to populations. The second panel from the top is the total signal at T = 0 and the 
remaining two panels per column are for T = 240(as) and T = 480(as). This form of the 
signal allows one to disentangle valence-core {uxa) transitions from Raman {uca) 

transitions. 
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FIG. 7: (Color online) (a): Resonant linear hybrid energy signal for coherences at multiple 
times T (top) and comparison of populations with initial coherences (bottom), (b): 

Resonant linear hybrid photon number signal for coherences at multiple times T (top) and 
comparison of populations with initial coherences (bottom). For time-dependence of 
coherences, T advances in units of ~60 attoseconds as one goes up the vertical axis. 


signals considered differ. We demonstrated the various signals for a model consisting of a 
few valence and a few core states and calculated the frequency-dispersed, photon number, 
and energy change signals for both linear or quadratic held intensity scaling and broad or 
hybrid spectral held shapes. 

There are two mechanisms whereby the photon modes can change their occupation num¬ 
bers. A photon can either redistribute amongst the modes, swapping from one frequency to 
another as in Raman interactions or, alternatively, it can be absorbed or emitted by matter. 
The former process is governed by the real parts of the polarizabilities a (and /3) while the 
latter is governed by the imaginary parts. In the oh-resonant regime, only the redistributive 
mechanism is operative. In the resonant regime, both contribute and oscillate with diherent 
phases in delay time T (the redistributive being a Sine funcion and the absorbtive/emittive 
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FIG. 8: (Color online) (Resonant linear (left colnmn) and quadratic (right column) hybrid 
photon number (top) and energy change (bottom) signals. These signals show G 
resonances at the Raman transitions Uac allowing their separation from the valence-core 

transitions 


being a Cosine). 

The spectrum of the nonstationary state (the energy levels of the occupied states), can be 
obtained by taking the Fourier transform of any of the time-dependent signals, resonant or 
off-resonant. The relative magnitude of the various peaks in these signals varies depending 
on the coupling constants for each signal. In the simpler, linear signals, the relative peak 
magnitudes are directly indicative of the product of the coherence pac and the polarizabil¬ 
ity aac, while in the quadratic case, there is an additional a and a free summation over 
the intermediate state. Using broadband pulses, we had to sacrihce the time resolution to 
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obtain frequency resolution of the ojac transitions. This limitation was surpassed by using 
a shaped, hybrid pulse consisting of broad and narrowband components. The resultant 
frequency-dispersed signals naturally show the peaks at the Uac transitions which oscillate 
in time. In the linear case, each peak oscillates with its own phase UacT + 0^^, thus allow¬ 
ing determination of the phase induced by the state preparation process. This technique 
therefore exploits the simultaneous resolution of narrowband and detected frequencies com¬ 
bined with the time-resolution granted by the broadband pulse to obtain the phases of the 
electronic wavepacket. Without frequency-dispersed detection, one is left with the photon 
number and energy change signals. Even though they contain less information than the 
frequency dispersed hybrid signals, the resolution granted by the narrowband pulse allows 
one to access essentially the same information as the frequency dispersed broadband signal. 

Finally, we note that the present formalism also applies to optical signals. X-ray signals 
allow one to probe valence excitations rather than vibrations. Furthermore, both the initial 
preparation and/or the preparation process can then be the product of an X-ray scattering, 
photoionization, or Auger process (as recently discussed in [24, 40]) in addition to an X-ray 
Raman process. The signals obtained by a Raman probe then detect both the amplitude and 
phase of the coherent superpositions of singly or doubly ionized states. These are signatures 
of many- body effects in the photoionization and the Auger processes. 

Compared to our previous work, here we present the most general systematic description 
of series of X-ray Raman signals. Using abbreviations of the signals of the present paper in 
our earlier work [15] we presented S^qv 2 i partially signals. Using the present 

terminology, the signals studied in [25] are Sq^ and . Finally, in the most recent work 
[26], we investigated signal. 


Appendix A: Multidimensional Off-Resonant Signals Initiated by an Impulsive Ra¬ 
man Process 

1. Two-pulse Time-Domain Experiment 

We now turn to a more elaborate experiment that probes the same material quantity 
described by Eqs. (22) but involves two broadband pulses separated in time by delay T 2 . 
At time Tq the system is prepared in p^c- Following the delay Ti = ti — tq the first probe 
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pulse Si centered at ti interacts with the system and after another delay T 2 = t 2 — ti the 
second pulse 82 centered at t 2 yields the two-dimensional spectra as a fnnction of the two 
time delays. The corresponding diagrams are given in Fig. 2(b). The difference is in the 
detection as well as the fact that the signal is linear in intensity of Si as well as 82 - The 
freqnency dispersed transmission is a 3D signal 


qUd) 

^Q12 


(ca, Ti, T2) 



X [\S2{uj - OJad)\\Sl{uJl - UJdc)\cOs{uJadT2 + 0ac(7'l)) 
- 1 ^ 2 ( 0 ; - Udc)\\Sl{Ui - UJad)\cOs{UdcT2 + 0ac(Ti))], 


(Al) 


the nnmber photon signal = 0 vanishes, and the energy change signal reads 




doJi 


\S2{uj)\\Si{uj)\^\pac\oif}a^ 

a,c,d 


da 


X [huJda\S-2{u: - UJad)\\Sl{u:i - Udc)\cOiS,{UadT2 + 0ac(Ti)) 
f)hJdc\S2{u: - UdMSiiUl - Uad)\cOs{UdcT2 + (t)ac{Tl))]. 


(A2) 


This is a snm of two flnxes a ^ d and c ^ d. Note, that compared to the signal (23), the 
two-pnlse signal (A2) carries information abont the phase and dynamics of the system after 
interacting with the probe Si dnring the delay T 2 . This is not accessible by the QB signal. 


2. Multidimensional Stimulated Raman Spectroscopy with an OfF-Resonant Ac¬ 
tinic Excitation 

There are multiple ways to prepare initial state pac snperposition. These include pho¬ 
toionization, Anger process, off-resonant or resonant impnlsive Raman excitation. So far 
we did not specify the preparation. In the standard formnlation of mnltidimensional spec¬ 
troscopy, the molecule is prepared by a series of short temporally well-separated pulses Sj 
at time r^-, j = 1,2,.... Followed by linear transmission of the probe (assnming no dynam¬ 
ics between preparation time Tq and interaction with probe at t) one can then expand the 



correlation function in Eqs. (2) - (5) 

((/|a®(t)e-s/^^-W'^^|p.(ro))) = \So{To)\^{[a{t),a{To)\, Pg) 

- ^|^o(T-o)n^(ri)n[a(t),a(ro)|,a(ri)],Pg) 

+ |^o(ro)ri^(ri)n^(r2)n[[a(t),a(ro)|,a(ri)],a(r2)],Pg) + ... (A3) 

The same result can be obtained by neglecting the preparation (at tq system is in pg) and 
expanding the dynamics between preparation time Tq and interaction with the probe pulse 
at t assuming series of short pulses prior to the probe pulse. 

In many X-ray spectroscopy applications, a single resonant actinic pulse £a prepares 
the system in the superposition of valence states. We consider an off-resonant preparation. 
We do not expand the signal in sum over states and obtain a general expression for the 
initial state described by the density matrix Pi(ro) assuming that the preparation process is 
temporally well separated from the probe 

_• fOO 

lft(ro))> = y y - r„)pQy(('')|p,)>. (A4) 

The linear frequency-dispersed transmission signals (7) then reads 

_07 ^OO 

^ J —oo J —oo 

(AS) 


Sl^\u,to,To) = 1 ^ I dt I d£\S{t-to)\^\SA{t'-To)\‘^{[a^^\t),a^^\t')],Pg), (A6) 


- 2 i 


r*oo /*oo 


h? 


' —oo o —oo 


r*oo /‘OO 


S'P(a;,to,To) =-X| / dt j dt'S*{t - to)S{t - to)\SA{t' - To)\‘^{[a^^\t),a^^\t')], pg), 


' —OO ^ —OO 


(A7) 

The quadratic frequency-dispersed transmission signal (17) is 

_o poo pt poo 

S^f^'^\uj,tQ,To) =X— I dt I dt' I dtV‘^(*"*°)T*(a;)T(t - to)|X(t'- to)nXA(^" - T-o)P 


OO J—OO J—OO 

X {[[a^^'>{t),a^^'>{t')],a^^'>{t")]pg). 


(AS) 




_o poo pt poo 

= / dt dt' / dt"\£{t - to)\^\£{t' - to)\^\£Ait" 

' ^ J —OO J —oo J —oo 


ro)P 


(A9) 
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O' poo pt poo 

slf\to,To) = dt dt' dt"S*{t-to)S{t-to)\S{t'-to)\^\SAit"-To)\‘^ 

^ J-oo J-oo J-oo 

x([[a(o)(t),a(°)(t')],a(°Ht")]p,). (AlO) 

These agree with original expressions obtained for homodyne detection [13]. The photon 
number and energy change signals can be obtained similarly. They provide a different gating 
of the same response function given by multiple commutators of the bare polarizability 

Appendix B: Off-Resonant Limit of the Resonant Polarizability 

We now examine the frequency-dispersed signal (31) and show that it reduces to the 
off-resonant result of Eq. (10). We first note that in the off-resonant case (Eqs. (32)) we 
can dehne the following quantity: 

= S*{uj)S{uj + uJac)al^J, (Bl) 

where 

^( 0 ) ■ ^a.x){G£ ■ Vxc) 

= 2^-ZT- ^^2) 

^xa 

X 

is the off-resonant polarizability. Choosing real basis set wave functions, the dipole moments 
are real. In this case the off-resonant polarizability (B2) is real as well. Note, that techni¬ 
cally this polarizability still contains field dependence via polarization of the probe pulse. 
Similarly 

= ^*(w)^(w - uJac)a^^- (B3) 

Substituting Eqs. (Bl) - (B3) into Eq. (31) we obtain (10). The other detection protocols 
and pulse shapes can be calculated similarly. 


Appendix C: Multidimensional stimulated Raman spectroscopy with resonant ac¬ 
tinic excitation 


If the preparation involves resonant actinic pulse the relevant excitation process is given 
by two diagrams shown in Fig. 9 and the straightforward calculation yields 


n = —aAA" 
Hac ^2 ^ac 


(Cl) 


30 



FIG. 9: (Color online) Set of ladder diagrams for resonant actinic preparation 


We can now substitute this matrix element into the signals and obtain a compact form in 
terms of a only. For the case of broadband pulse Eqs. (31), (33) and (35) yield 



) . (C2) 



(C3) 
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